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■ ABSTRACT 

. I A large fraction of brown dwarfs and low-mass hydrogen-burning stars may form 

■ by gravitational fragmentation of protostellar discs. We explore the conditions for 

disc fragmentation and we find that they are satisfied when a disc is large enough 

Q , (> 100 AU) so that its outer regions can cool efficiently, and it has enough mass 

5— ( ' to be gravitationally unstable, at such radii. We perform radiative hydrodynamic 

[/3 I simulations and show that even a disc with mass 0.25 Mq and size 100 AU fragments. 

. The disc mass, radius, and the ratio of disc-to-star mass (Md/M^ « 0.36) are smaller 

than in previous studies (Stamatellos & Whitworth 2009a). Wc find that fragmenting 
I discs drastically decrease in mass and size within a few 10^ yr of their formation, since 

^ I a fraction of their mass, especially outside ^ 100 AU is consumed by the new stars and 

ly^ ' brown dwarfs that form. Fragmenting discs end up with masses ~ 0.001 — 0.1 Mq, and 

ly^ , sizes ^ 20 — 100 AU. On the other hand, discs that are marginally stable evolve on 

■rj" ■ a viscous timescale, thus living longer (~ 1 — 10 Myr). We produce simulated images 

CO ' of fragmenting discs and find that observing discs that are undergoing fragmentation 

. is possible using current (e.g. IRAM-PdBI) and future (e.g. ALMA) interferometers, 

but highly improbable due to the short duration of this process. Comparison with 
observations shows that many observed discs may be remnants of discs that have 
fragmented at an earlier stage. However, there are only a few candidates that are 
possibly massive and large enough to currently be gravitationally unstable. The rarity 
of massive (> 0.2 Mq), extended (> 100 AU) discs indicates either that such discs 

■ are highly transient (i.e. form, increase in mass becoming gravitationally unstable 
?-H ' due to infall of material from the surrounding envelope, and quickly fragment), or 
^ that their formation is suppressed (e.g. by magnetic fields). We conclude that current 

observations of early-stage discs cannot exclude the mechanism of disc fragmentation. 

Key words: Stars: formation - Stars: low-mass, brown dwarfs - accretion, accretion 
disks - Methods: Numerical, Radiative transfer. Hydrodynamics 



1 INTRODUCTION 

The formation of low-mass hydrogen-burning stars, brown 
dwarfs (BDs) and planetary-mass objects is a critical prob- 
lem, since it addresses the question of why there is a mini- 
mum mass for star formation. These objects are difficult to 
form by gravitational fragmentation of unstable g for 
masses in the BD regime (J; 80 Mj, where Mj is the mass of 
Jupiter) a high density is required for the gas to be Jeans un- 
stable. Thus, if BDs form according to the standard model of 
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star formation, i.e. the collapse of a prestellar core, then the 
pre-BD core must have a density of 10~^®g cm~'^. Padoan 
& Nordlund (2004) and Hennebelle & Chabrier (2008) sug- 
gest that these high density cores may form by colliding flows 
in a turbulent magnetic medium. However, it has not been 
demonstrated that this mechanism works when the thermo- 
dynamics and hydrodynamics are properly simulated, and 
there is no obvious reason why it should reproduce the ob- 
served binary properties of BDs and low-mass stars. 

An alternative way to reach the high densities required 
for the formation of BDs is in protostellar discs (Whitworth 
& Stamatellos 2006; Stamatellos et al. 2007b). These discs 
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form around newly born stars and grow quickly in mass by 
accreting material from the infalling envelope (Terebey et al. 
1984; j0rgensen et al. 2009; Vorobyov & Basu 2010; Machida 
& Matsumoto 2010; Walch et al. 2010). They become un- 
stable and fragment if the matter infalling on them cannot 
efficiently redistribute its angular momentum outwards in 
order to accrete onto the central star (Attwood et al. 2009). 
Whitworth & Stamatellos (2006) argue that fragmentation 
happens in the outer disc where the cooling time is short 
enough (Matzncr & Levin 2006; Bolcy 2009). Stamatellos 
et al. (2007b) show that this is a viable mechanism, and 
Stamatellos & Whitworth (2009a) determine the statistics 
of the low-mass objects produced by this mechanism, per- 
forming an ensemble of radiative hydrodynamic simulations 
of fragmenting discs. The disc fragmentation model repro- 
duces the critical constraints set by the observed statistical 
properties of low-mass objects, and these have not yet been 
explained by other formation mechanisms. More specifically 
the model reproduces the shape of the mass distribution of 
low-mass objects (IMF), the lack of BDs as close compan- 
ions to Sun- like stars (i.e. the BD desert), the presence of 
discs around BDs, the statistics of low-mass binary systems, 
and the formation of free-floating planetary mass objects 
(Stamatellos & Whitworth 2009a). 

The main argument against the formation of low-mass 
stars and BDs by fragmentation of discs is whether the mas- 
sive (a few times 0.1 Mq), extended discs (^ 100 AU) as- 
sumed by this mechanism can be realized in nature. From 
a theoretical point of view the formation of such discs is in- 
evitable considering the angular momentum content of their 
parental cores. For example, a 1.4 Mq prestellar core with ra- 
tio of rotational to gravitational energy /3 = TZ/\il\ will - if it 
collapses monolithically - form a protostellar disc with outer 
radius R^.^c ~ 400 AU (/3/0.01). The observations of Good- 
man et al. (1993) indicate that many prestellar cores have 
/3 0.02, hence the formation of extended, self-gravitating 
discs should be rather common. 

Recent numerical studies reach mixed conclusions on 
this topic. Simulations that ignore the effects of magnetic 
fields (Attwood et al. 2009; Machida & Matsumoto 2010) 
find that such large discs form frequently in simulations of 
turbulent and/or rotating molecular clouds. The presence 
of magnetic fields is expected to act against the formation 
of a centrifugally supported disc because angular momen- 
tum is removed by magnetic effects (e.g. magnetic braking, 
outflows). However, it is uncertain whether, in realistic con- 
ditions, magnetic fields can totally suppress the formation 
of self-gravitating discs. Hennebelle & Promang (2008) and 
Henncbelle & Tcyssicr (2008) find that in the ideal MHD 
approximation the formation of a disc is suppressed if the 
magnetic field is strong enough (with strength similar to 
what is inferred by observations, e.g. Crutcher et al. 2004, 
Falgarone et al. 2008) and parallel to the rotation axis of 
the collapsing star-forming core. If the magnetic field is mis- 
aligned then its effect is weakened and discs may form, but 
fragmentation seems to be suppressed (Hennebelle & Cia- 
rdi 2009). These results axe supported by ideal MHD sim- 
ulations that include the effects of radiative transfer (Com- 
mcrcon et al. 2010). Disc formation and fragmentation in 
strongly magnetized clouds appears to be possible only for 
rapidly rotating clouds (Machida et al. 2005). Ambipolar 
diffusion tends to weaken the effects of the magnetic fields. 



but for realistic magnetic field strengths and ionization rates 
not enough to limit the effect of magnetic breaking (Mellon 
k Li 2008, 2009). The situation changes in resistive MHD 
calculations. Ohmic dissipation is probably important at 
high densities and may dominate over other processes at 
densities 10^^ cm"^ (Nakano et al. 2002). Machida et al. 
(2008, 2010) using a resistive MHD method, for a wide range 
of rotational-to-magnetic energies, find that disc formation 
is possible. In these simulations fragmentation happens ei- 
ther on a large scale and produces wide binaries, or after 
the formation of the first core and produces close binaries. 
Krasnopolsky et al. (2010) find that a resistivity on the or- 
der of 10^" cm"^ s^^ is needed for the formation of discs with 
sizes on the order of 100 AU. 

The goal of this paper is to determine the lower limits 
of fragmentation in terms of disc mass and size, and com- 
pare these with observations of early-stage discs. Two condi- 
tions have to be satisfied for a disc to fragment: (i) the disc 
must be Toomrc unstable i.e. it has to be massive enough so 
that gravity can overcome thermal and local centrifugal sup- 
port (Toomre 1964) , and (ii) the disc must cool efficiently so 
that proto-condensations forming in the disc do not simply 
undergo an adiabatic bounce and dissolve (Gammie 2001). 
These criteria can be satisfied for discs considerably smaller 
in mass and size than the discs assumed by Stamatellos & 
Whitworth (2009a). In particular, the fragmentation crite- 
ria are fulfilled for discs that (i) are large enough so that 
their outer regions can cool fast enough, i.e. for discs with 
radii ^ 100 AU, and (ii) have enough mass (i.e. an appro- 
priate surface density profile) so that they are Toomre un- 
stable at such radii. In this paper we present an ensemble 
of simulations exploring the possibility of fragmentation of 
small, marginally unstable discs. We find that even discs 
with masses 0.25 Mq and sizes 100 AU fragment. We pro- 
duce simulated images of such discs that show that they are 
observable by current facilities (e.g the IRAM Plateau de 
Bure interferometer; hereafter PdBI), although they should 
be rarely observed as they quickly fragment. Gomparison 
with observations shows that a few marginally unstable discs 
may exist. Additionally, many discs are consistent with be- 
ing remnants of disc that have fragmented in the past. 

The paper is structured as follows. In Section 2 we de- 
scribe the initial conditions of the disc simulations, the com- 
putational method used, and the results regarding the lower 
limits of disc fragmentation as they are derived by the sim- 
ulations. In Section 3 we produce synthetic images of discs 
during the fragmentation process, and in Section 4 we com- 
pare the properties of fragmenting discs with observations 
of candidate early-stage discs. Finally, in Section 5 we sum- 
marise this work and discuss its implications. 



2 HYDRODYNAMIC SIMULATIONS OF 
MARGINALLY UNSTABLE DISCS 

2.1 Disc initial conditions: Density and 
temperature profiles 

We assume a star-disc system in which the central primary 
star has initial mass = 0.7 Mq. The initial disc mass, 
Mp, and the initial disc radius, JJ^, axe free parameters. 
The disc mass is varied from 0.15 to 0.30 Mq and the disc 
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radius from 60 to 150 AU (see Table [1}. The initial surface 
density is 

E„(i?) = E(lAU) (A)'" , (1) 
and the disc temperature 

T„(i?) = 250 K (A) VlOK, (2) 

where E(IAU) is determined by the disc mass and radius, p 
and q are free parameters, and R the distance from the cen- 
tral star measured on the disc midplane. The values of q and 
p are rather uncertain and may evolve with time. Andrews et 
al. (2009) observed a sample of circumstellar discs in Ophi- 
uchus and estimated that the disc surface density drops with 
an exponent p ~ 0.4—1.0. In Taurus- Auriga and Ophiuchus- 
Scorpius star formation regions Andrews & Williams (2007) 
find a median p ~ 0.5; they also find a temperature profile 
with g Ri 0.4 - 0.74. Osterloh & Beckwith (1995) find that 
the disc temperature drops with an exponent q « 0.35 — 0.8. 

Previous studies (e.g. Stamatellos & Whitworth 2009a) 
have used p = 1.75 and q = 0.5, resulting in a disc with 
almost uniform Toomre parameter Q. The conditions are 
more favourable for fragmentation for lower values of p (as 
the disc has more mass at larger radii), and higher values of 
q (as the temperature drops faster with the distance from 
the central star). In this study we set p — 1, hence there 
is equal mass per unit radius in the disc, and q — 0.75 or 
q — 0.5 (see Table[T|. These parameters make fragmentation 
easier than in previous simulations (Stamatellos et al. 2007; 
Stamatellos & Whitworth 2009a), but they are not extreme 
values when compared with disc observations. 

2.2 Computational method 

We use the SPH code seren (Hubber et al. 2011), which in- 
vokes an octal tree (to compute gravity and find neighbours) , 
adaptive smoothing lengths, multiple particle timesteps, 
and a second-order Leapfrog integration scheme. The code 
uses time-dependent artificial viscosity (Morris & Monaghan 
1996) with parameters a* — 0.1, /? — 2a and a Balsara 
switch (Balsara 1995), so as to reduce artificial shear viscos- 
ity. 

The radiative processes that regulate the disc thermo- 
dynamics affecting the fragmentation process are treated 
with the Stamatellos et al. (2007a) diffusion method that 
includes radiation heating from the central protostar as de- 
scribed in Stamatellos et al. (2009b) and in Forgan et al. 
(2009). The method takes into account compressional heat- 
ing, viscous heating, radiative heating by the background, 
and radiative cooling. It performs well, in both the optically 
thin and optically thick regimes, and has been extensively 
tested (Stamatellos et al. 2007a; Forgan et al. 2009). In par- 
ticular it reproduces the detailed 3D results of Masunaga & 
Inutsuka (2000), Boss & Bodenheimer (1979), Boss & My- 
hill (1992), Whitehouse & Bate (2006), and also the analytic 
resuhs of Spiegel (1957) and Hubeny (1990). 

The gas is assumed to be a mixture of hydrogen and 
helium. We use an EOS (Black & Bodenheimer 1975) that 
accounts (i) for the rotational and vibrational degrees of free- 
dom of molecular hydrogen, and (ii) for the different chem- 
ical states of hydrogen and helium. We assume that ortho- 



Table 1. Parameters of the simulations of marginally unstable 
discs. MJ^: disc initial mass, R^ : disc initial radius, M^: disc final 
mass (i.e. at f = 40 kyr), R^: disc final radius, p: exponent of the 
disc surface density profile, q: exponent of the disc temperature 
profile, A'^f; number of objects formed. The colour corresponds to 
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Table 2. The properties of the objects formed by disc fragmen- 
tation. We list their final mass (Mf), formation radius (rjn), final 
status (i.e. whether they remain bound to the central star or are 
ejected from the disc), and their final semi- major axis (afin) and 
eccentricity (e) (if they are still bound to the central star), or 
their velocities (foj) (if they are ejected from the disc). 
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and para-hydrogen are in equilibrium. For the dust and gas 
opacity we use the parameterization by Bell & Lin (1994), 
k{p,T) = Ko p" , where kq, a, b are constants that de- 
pend on the species and the physical processes contributing 
to the opacity at each p and T. The opacity changes due to 
ice mantle melting, the sublimation of dust, molecular and 
H" contributions, are all taken into account. 

Sinks are invoked in the simulations in order to avoid 
very small timesteps when the density becomes high. Sinks 
are created when a particle reaches a critical density of 
psink = 10~^g cm"' and the particles in its neighbourhood 
(i.e. its 50 neighbours) are bound. The sink radius is set to 

1 AU. Once a particle is within this radius and bound to the 
sink then it is accreted onto it. 

The disc is represented by 1.5 x 10"" SPH particles, 
which means that the minimum resolvable mass (corre- 
sponding to the number of neighbours used, i.e. 50 SPH 
particles) is f» 0.1 Mj for the most massive disc simulated 
(Md = 0.3 Mq). Meru & Bate (2010) find that disc frag- 
mentation simulations have probably not converged using 

2 X 10® particles. However, their study is based on a pa- 
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Figure 1. Toomre parameter (top) and cooling time in units of the local orbital period (bottom) for discs that do not fragment, at 
1000 yr from the beginning of the simulation. Both are azimuthally averaged. The different colours correspond to different runs as listed 
in Table [1] 
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Figure 2. Toomre parameter (top) and cooling time in units of the local orbital period (bottom) for discs that fragment, at 1000 yr from 
the beginning of the simulation. Both are azimuthally averaged. The different colours correspond to different runs as listed in Table [T] 
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rameterized cooling prescription, and on only two high res- 
olution (with 1.6 X 10^ particles) simulations; more work is 
needed to establish the effect of resolution in disc fragmen- 
tation. Stamatellos & Whitworth (2009b) have shown that 
1.5 X 10^ SPH particles are adequate to properly resolve the 
Jeans mass, the Toomre mass, and the vertical structure of 
the disc (Nelson 2006). Additionally, Stamatellos & Whit- 
worth (2008) demonstrate that, with 1.5 x 10^ particles, the 
prescription for the radiative transfer used in this work, re- 
produces well the vertical disc temperature profile derived 
analytically by Hubeny (1990). 

2.3 The lower limits of disc fragmentation 

We examine the possibility of fragmentation for discs with 
masses from 0.15 to 0.30 Mq and radii from 60 to 150 AU 
(around a 0.7 Mq star) (see Table [l}. All the discs are 
Toomre unstable or close to being so (Figs. [Tl [2]). Gravi- 
tational instabilities develop and grow in all simulations; in 
some of the simulations these instabilities become stronger 
leading to disc fragmentation within a few thousand years 
(Fig. m first column), and in other simulations they weaken 
resulting in an almost axisymmetric disc in a self-regulating 
state (Lodato & Rice 2004, 2005). The simulations are fol- 
lowed for long enough time (40 kyr or equivalently ~ 20 
orbital periods at 150 AU) to ensure that gravitational insta- 
bilities have adequate time to grow. We find that even a disc 
with mass of 0.25 Mq and radius of 100 AU fragments; these 
values and the ratio of disc-to-star mass {Mo /M* « 0.36) 
are considerably lower than the ones assumed by Stamatellos 
& Whitworth (2009a) (0.7 Mq, 400 AU, and Md/M* = 1 
respectively) . 

A common characteristic of all simulations is that the 
the gravitational instabilities provide an effective way to re- 
distribute angular momentum. As a result the discs spread 
from an initial size of 60-100 AU to more than 200 AU. The 
discs can cool efficiently only outside 70-100 AU (e.g. Sta- 
matellos & Whitworth 2009b). If the disc remains Toomre 
unstable as it spreads then it fragments (runs 1, 3, 5; Fig. (2]). 
For example, in run 5 (see Table [T]), the disc has an initial 
size of 100 AU, but within 1000 yr its size has reached 250 
AU. Subsequently, the disc fragments to form 4 bound ob- 
jects at distances 100, 140, 175, and 200 AU (see Table[21). A 
similar process of disc expansion and fragmentation is sug- 
gested by Clarke (2009) but in her work the disc spreading 
happens on a viscous timescale, i.e. it takes much longer (a 
few 10^ yr). 

Typically the outcome of disc fragmentation are multi- 
ple low-mass objects, i.e. low-mass stars, BDs and planetary- 
mass objects (see Table [2]). Similar to the Stamatellos & 
Whitworth (2009a) simulations, the object that forms first 
migrates inwards in the inner disc region and increases in 
mass by accreting material, to eventually become a low- 
mass hydrogen-burning star. Objects that form later, farther 
out in the disc, have lower mass, i.e they are usually BDs 
and planetary-mass objects. In the three cases of fragment- 
ing discs 1, 3 and 4 objects form in each disc, respectively. 
Where multiple objects form in the disc, they dynamically 
interact with each other and the central star, and some of 
them are ejected from the system. This is what happens to 
the only planetary-mass object formed in the simulations. 
These results are qualitatively similar to the ones of Sta- 
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Figure 3. Disc mass (top) and disc radius (bottom) for frag- 
menting (dashed lines) and non-fragmenting discs (solid lines) 
at 3 times (0 kyr: black/square, 15 kyr: blue/circle, 40 kyr: 
red/triangle). The fragmenting discs reduce in mass and size 
quickly, whereas the non-fragmenting discs expand and dissipate 
slowly (see text). 



matellos & Whitworth (2009a), but a larger ensemble of 
simulations is needed for a more quantitative comparison. 
The main difference between the two sets of simulations is 
that fewer objects are formed in smaller discs (i.e. 1 — 4), 
as compared with larger discs (i.e. 5 — 12, Stamatellos & 
Whitworth 2009a). 

Following fragmentation, the disc around the central 
star is considerably reduced in mass and size (Table [TJ 
Fig. 131), because a fraction of its mass is going into the new 
objects that form in its outer regions. The three fragment- 
ing discs (runs 1, 3, 5) end up with mass 0.002 — 0.07 Mq, 
and size 20 - 70 AU, after only 40 kyr (Table [T] Fig. O. 
On the other hand, in discs that do not fragment the re- 
distribution of angular momentum results in the inner disc 
regions accreting onto the star and the outer regions expand- 
ing (Lynden-Bell & Pringle 1974). Initially this happens fast 
since the gravitational torques provide a sufficient effective 
viscosity, but then it slows down as the disc spreads and 
becomes more stable. This process is probably not affected 
by the presence of magnetic fields as the magneto-rotational 
instability that can provide an additional source of viscosity 
is presumed to be activated only in the inner disc regions 
(within a few AU from the central star; e.g. Zhu et al. 2009). 
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Figure 4. Hydrodynamic simulations (first column) and simulated 1.3-mm observations of fragmenting discs at the distance of Taurus 
(140 pc; second column), and Perseus (250 pc; third column). 6 snapshots are plotted. The simulations depict the logarithm of column 
density (g cm~^). The synthetic images refer to flux density (mjy/beam). The dashed contours on the synthetic images show the typical 
5(T detection level (~3 mJy/beam) that can be achieved by the PdBI (A conflguration, i.e. PdB-A). The white contours are levels of 20(t, 
and the following black contours correspond to 50, 100 and 200(7 levels. The ellipse in the bottom-left corner of each simulated image 
represents the FWHM beam size (resolution element) obtained with PdB-A. 
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2.4 Even lower limits for fragmentation? 

The simulations presented here explore only a small fraction 
of the parameter space. The criteria for disc fragmentation 
could be satisfied for even smaller disc masses and sizes. 
More specifically, the limits for fragmentation could be lower 
in the following cases. 

(i) For a less steep disc density profile. We have assumed 
that p = 1, but observations of T Tauri discs suggest that p 
could be as low as ~0.4 (see Section [2.11) . This means that 
there could be more mass at larger disc radii. The value of 
p for early-stage discs is uncertain. 

(ii) For a steeper disc temperature profile. We have as- 
sumed (in most simulations) that q = 0.75, whereas obser- 
vations of T Tauri discs suggest that q could be as high as 
~ 0.8 (see Section 12. ip . This results in cooler, and hence 
more unstable discs. The value of q for early-stage discs is 
also uncertain. 

(iii) For a lower mass central protostar. We have assumed 
= 0.7 Mq, but at the initial stages of their formation 

protostars are less massive. This has two consequences: (a) 
their discs could be unstable even if they have masses below 
~0.01 Mq (e.g. Inutsuka et al. 2010), and (b) their discs 
can cool fast enough even if they are small, as the minimum 
radius for disc fragmentation scales as A/^^ (Whitworth & 
Stamatellos 2006; Clarke 2009). 

(iv) For discs that are perturbed by passing by stars. It 
has been shown that such encounters promote disc fragmen- 
tation (Thies et al. 2010). 

Thus, the limits (mass, size) of disc fragmentation are 
most probably lower than the ones reported here. In a future 
publication we will estimate the disc fragmentation limits 
considering the above factors. 



3 SYNTHETIC IMAGES OF FRAGMENTING 
DISCS 

Comparing the results of current interferometric millimetre 
observations of young protostars (Class 0, Class I objects) 
with the predictions (e.g. mass and size of the discs, multi- 
plicity of the stars formed) of the simulations is important 
for understanding the first steps of low-mass star formation 
by disc fragmentation. In the following paragraph we de- 
scribe the method used to produce synthetic images from 
numerical simulations of fragmenting discs, so as to deter- 
mine whether such discs and their spiral structure are de- 
tectable by current facilities, and to compare these discs with 
observations. 

The column density maps produced by the numerical 
simulations are scaled to the distance of the Taurus (140 pc) 
and Perseus (250 pc) star forming regions. Then, they are 
transformed to flux density images (in mjy/beam) assum- 
ing typical dust properties for Class objects: optically thin 
emission at 1.3 mm, a dust temperature Tdust = 10 K and a 
dust emissivity Ki.smm = O.Olcm^g"^. (Note that the dust 
temperatures may be higher and that the dust opacity is un- 
certain by a factor of at least ~ 2). Finally, the resulting im- 
ages are adapted to the GILDAqjdata format, in order to be 



processed with the PdBI simulator, mostly comprising the 
'uvjnodel' and 'uvjnap' tasks. These tasks allow us to pro- 
duce simulated 1.3 mm continuum images, i.e. images that 
would be obtained if these objects were observed with PdBI. 
During this last step, we first create the synthetic reference 
uv-coverages corresponding to the PdBI in its most extended 
(A) configuration. Pointing towards the coordinates of the 
source is also considered; then 'uvjnodeP convolves the flux 
density maps with the reference uv-coverages, and finally 
'uv_map' inverts the resulting uv-tables to create the syn- 
thetic 1.3 mm continuum maps. 

Fig. [4] shows the synthetic images for 6 snapshots of 
simulation 3 (see Table [l]) (second and third column) to- 
gether with the logarithmic column density maps from the 
simulation (first column). The images are 1.3 mm contin- 
uum images corresponding to simulated observations us- 
ing the A-array of the PdB interferometer at Taurus (sec- 
ond column) and Perseus (third column) . The dashed white 
contours in the image correspond to the typical 5a level 
3 mJy/beam) that can be achieved in the 230 GHz 
continuum with the A-array of PdBI, taking into account 
the increased sensitivity of the newly installed WideX cor- 
relator. We see that the PdBI is sensitive enough to detect 
fragmenting discs as extended structures (expected angular 
diameter ~ 2"-3" compared to a PdBI beam size ~ 0.5" 
in the A configuration). On the other hand, PdBI does not 
provide enough resolution to discern any spiral structure in 
the disc. However, the upcoming ALMA should in princi- 
ple be able to detect such spiral structures out to distances 
of the Taurus- Auriga star-forming regions, i.e. at ~ 140 pc 
(Cossins et al. 2010). 

The disc around the protostar (Fig. |4]) initially is large 
(100 AU) and massive (0.3 Mq) but within a few thousand 
years it reduces in size (^ 35 AU) and mass (~ 0.03 Mq), 
as it fragments to form companions (i.e. low-mass stars and 
BDs). Thus, though fragmenting discs can in principle be 
detected, they may be rarely observed due to their short 
lifespan (see next Section). Companions to the central proto- 
star (which are signposts of fragmentation) can be detected 
if they are attended by discs that have enough mass. When 
such discs are disrupted (e.g. by dynamical interactions with 
other companions in the disc) then the low-mass companions 
may remain undetected (e.g. 4"^ row, 3"^^^ column) . They also 
have to be away from the central protostar, so that they can 
be resolved as separate entities (e.g. in 5**^ row, 2'^'* column, 
the close companion is blended with the central protostar). 



4 COMPARISON WITH DISC OBSERVATIONS 

4.1 Uncertainties in calculating disc masses 

The masses and radii of early stage discs are difficult to de- 
termine. They are rather uncertain, maybe up to an order 
of magnitude (see discussion in j0rgensen et al. 2009). Discs 
around newly formed protostars (i.e. Class 0) are generally 
hidden within the envelope that surrounds the protostar, so 
that radiative transfer models are needed in order to dis- 
tinguish the submm and mm emission of the disc from the 



^ Grenoble Image and Line Data Analysis Sys- 



tem, software provided and developed by IRAM 
( http :77m . Iram . f r/IRAHFR/GILDAS | 
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emission of the envelope. For example, the j0rgensen et al. 
(2009) models suggest that in Class objects 68% of the 
observed interferometric flux is from the disc (87% in Class 
I objects). Once the integrated flux, Fint, from the disc is 
known then its mass is computed using 

Md = 4.0 X 10^ Mq (^77^ V X 

^ \^ 10 mjy J \ 140 pc J 

0.01 cm2 g-i ) V 20K J • ^'^> 

Typical values for the dust temperature and the opacity are 
< Tdust >= 20 K and Ki.smm = 0.01 cm^ g~^, respectively 
but they are rather uncertain. Ossenkopf & Henning (1994) 
estimate that opacities may increase as much as an order of 
magnitude depending on the density. The dust temperatures 
in discs are also poorly constrained and the use of an average 
dust temperature may be inappropriate. 



4.2 Candidate unstable discs 

The fragmenting discs (runs 1, 3, 5; Table [T]) in our sim- 
ulations reduce to a mass of 0.002 — 0.07 Mq, and a size 
of 20 - 70 AU, within 40 kyr (see Fig. [3] and Table [l]). 
This happens so fast, that large, unstable discs or discs that 
are currently undergoing fragmentation should be rarely ob- 
served (see Discussion). There are only a few observations 
of discs that (considering the uncertainties in the estimated 
disc masses and radii, see Section [4. 11 and in disc fragmen- 
tation models, see Section [2. 4|) . may be gravitationally un- 
stable. A recent example is the disc around HL Tan, a Class 
II object, which has a mass of around 0.1 Mq and extends 
to at least 100 AU, with a possible low-mass fragment at 
~ 65 AU (Greaves et a. 2008). It has been shown that this 
fragment may have formed by disc fragmentation (Greaves 
et a. 2008). Another example of a possible unstable disc, 
is Serpens FIRS 1 with mass of 1.0 Mq and radius 300 AU 
(Enoch et al. 2009a); this however may be a rotating collaps- 
ing spheroid. Other discs that may be gravitationally unsta- 
ble are GSS39 (Class II; 0.143 Mq, 200 AU) and DoAr 25 
(Class II: 0.136 Mq, 80 AU) in Ophiuchus (Andrews et al. 
2009), and HC189 (0.25 Mq, < 100 AU) in the Orion Neb- 
ula Cluster (Eisner et al. 2008). Another candidate is the 
disc around IRAS 16293-2422B (Class 0; Rodriguez et al. 
2005) which is resolved with the VLA and has an estimated 
mass of 0.3 - 0.4 Mq and size of 26 AU. 



4.3 Remnants of discs that may have fragmented 

Discs that have masses and radii consistent with being rem- 
nants of discs that have fragmented in the past, are rou- 
tinely observed at the Class I and Class II stages; how- 
ever, it is difficult to establish whether fragmentation has 
occurred. Discs with masses and sizes consistent with the 
final {t — 40 kyr) masses and sizes in our simulations 
(0.002 - 0.07 Mq, 20 - 70 AU) have been observed (e.g. 
Eisner et al. 2005, 2006; 2008; Rodriguez et al. 2005; Eisner 
& Carpenter 2006; Brinch et al. 2007; Andrews & Williams 
2007; Andrews et al. 2009; J0rgensen et al. 2009; Brinch et 
al. 2009; Tobin et al. 2010). 



4.4 Surveys of early-stage discs 

Observations of early-stage discs are difficult as these discs 
are deeply embedded in their parental clouds. In the follow- 
ing paragraphs we discuss a few surveys that report rela- 
tively massive disc candidates. 

j0rgensen et al. (2009) using submm/mm observations 
and radiative transfer modelling computed disc masses from 

0. 017 to 0.46 Mq, and from 0.008 to 0.053 Mq in a sam- 
ple of Class and Class I sources, respectively. All of these 
discs were unresolved or marginally resolved (size J: 1" — 3", 

1. e. ~ 125 — 700 AU). In a few of the Class I sources of 
the j0rgensen et al. sample HCO"'" 3 — 2 line observations 
revealed signs of Keplerian rotation confirming their disc 
status (e.g. L1489-IRS; Brinch et al. 2007). Scattered light 
images of L1527 and radiative transfer modelling infer a disc 
with radius ~ 190 AU and mass of 0.005 Mq (Tobin et al. 
2010). We note however that the mass of this disc could be 
up to 0.1 Mq if different opacities are assumed (Loinard et 
al. 2002). Interestingly, a close companion at projected dis- 
tance of ~ 24 AU is also present (Loinard et al. 2002), indica- 
tive of disc fragmentation. In the same sample of J0rgensen 
et al. (2009), NGC1333-IRAS2 (Class 0) is estimated to have 
a radius of _R = 150 AU and mass of M = 0.1 Mq although 
no sighs of Keplerian rotation were found (j0rgensen et al. 
2005; Brinch et al. 2009). The estimated masses and sizes 
L1527 and NGC1333-IRAS2 are contested as recent 1.3-mm 
observations (Maury et al. 2010) reveal no discs larger than 
the resolution limits of PdB-A (-- 100 AU). 

Eisner et al. (2005, 2006, 2008) observed young proto- 
stars (Class Os and Class Is) in the Orion Nebula Cluster and 
using radiative transfer models inferred discs with masses up 
to 0.39 Mq and radu up to 220 AU. 

Andrews et al. (2009) used the SMA to achieve high 
angular resolution of ^ 0.3", i.e. ~ 40 AU in Ophiuchus and 
observed discs with masses up to 0.14 M© sizes up to 200 
AU; some of them show signs of Keplerian rotation (evident 
in CO j = 3 — 2 line observations) . In the Taurus- Auriga 
and Ophiuchus-Scorpius star forming regions Andrews & 
Williams (2007) found disc masses of up to 0.17 Mq and 
sizes up to 700 AU (most of these around Class II objects) . 

4.5 Remarks 

Discs are probably present at the earliest phases of star for- 
mation, i.e. in Class and Class I objects. However, their 
masses and sizes are debated. There are a few candidates 
that they may be relatively massive (a few 0.1 Mq) and 
large 100 AU), but current observations are not conclu- 
sive (e.g. see Maury et al. 2010). Comparing typical masses 
and radii derived from disc observations with those of frag- 
menting discs in Table [1] (see also Fig. |3} we find that most 
of these discs are not gravitationally unstable. However, it 
cannot be excluded that these discs were gravitationally un- 
stable and that have fragmented in the past (i.e. in last 
~ 40 kyr). 



5 DISCUSSION 

The ensemble of disc simulations presented here show that 
protostellar discs fragment when they are large enough 
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100 AU) so that their outer regions can cool efSciently, 
and massive enough to be gravitationally unstable, at such 
radii. We find that even a disc with moderate size (100 AU) 
and moderate mass (0.25 Mq) around a 0.7-Mq protostar 
(Md/M* ~ 0.36) fragments. These values are considerably 
smaller than those in previous studies; Stamatellos & Whit- 
worth (2009a) simulated discs with size 400 AU and mass 
0.7 Mq, around a 0.7-Mq protostar (Md/M* = 1). Even 
less massive and smaller discs may be able to fragment (see 
Section [2. 4p : a larger parameter space needs to be explored 
in order to determine the lowest limits of disc fragmentation. 

The properties of the objects formed by the fragmen- 
tation of the discs presented in this paper show similarities 
with the properties of the objects produced in larger, more 
massive discs (Stamatellos & Whitworth 2009a). Most of the 
objects formed are BDs, but low-mass H-burning stars, and 
planetary-mass objects also form. Low-mass H-burning stars 
generally orbit closer to the central star, BDs are on wider 
orbits (170—1250 AU), whereas planetary-mass objects tend 
to be ejected from the disc, due to dynamical interactions. 
The most notable difference is that fewer objects form in 
smaller discs. A larger ensemble of simulations is needed to 
make a more quantitative comparison with the results of 
Stamatellos & Whitworth (2009a). 

Fragmenting discs can in principle be detected by cur- 
rent (e.g. PdBI; Section[3]) and future (ALMA; Cossins et al. 
2010) interferometers. However, the chance to observe such 
discs during fragmentation is small, as fragmentation hap- 
pens in a dynamical timescale, i.e. within a few thousand 
years. Hence, a large number of protostars need to be ob- 
served. Based on first results (Andre et al. 2010), the total 
number of Class protostars that are expected to be found 
in nearby clouds by the Herschel Gould Belt survey is 350. 
Assuming that (i) all protostars are attended by discs that 
grow, become unstable, and fragment during the Class 
phase, (ii) the lifetime of the Class phase is ~ 1.5 x 10^ 
yr (Evans et al. 2009; Enoch et al. 2009b), and (iu) in the 
disc fragmentation scenario the disc fragments and therefore 
dissipates within 1.5 x 10* yr (see Fig. O red triangles), it 
follows that ~ 10% of Class objects could have extended 
unstable discs. However, Stamatellos & Whitworth (2009a) 
find that typically a few BDs form in each fragmenting disc 
around a Sun-like star, and suggest that only ~ 20% of Sun- 
like stars need to have unstable extended discs to produce 
all BDs and a large fraction of low-mass hydrogen-burning 
stars. Hence, assuming that disc fragmentation delivers half 
of the observed BDs, i.e. that only ~ 10% of Sun-like stars 
have discs that fragment, then the percentage of Class 
objects with unstable extended discs is only ~ 1%. This per- 
centage could be even lower if a smaller fraction of BDs 
forms by fragmentation of early stage discs or if discs need 
not be that extended to fragment (see Section [2. 4[) . On the 
other hand it could be higher if the lifetime of the Class 
phase is shorter and the disc dissipation timescale longer. 
For example, assuming a Class lifetime of ~ 3 x 10* yr 
(Andre et al. 2000) and a similar disc dissipation timescale 
(•^ 3 X 10* yr) then the percentage of Class objects with 
large extended discs is ~ 10% (for half of the observed BDs 
to form by disc fragmentation). These estimates may also 
depend on how fast a disc assembles around a young pro- 
tostar forming in a collapsing molecular cloud core; larger 
scale simulations (Walch et al. in prep.), including the effects 



of radiative transfer and magnetic fields (e.g. CommerQon et 
al. 2010), are needed for assessing this. 

It is possible to observe discs that are remnants of disc 
that have fragmented in the past. Our simulations show that 
such discs may reduce to sizes ~ 20 — 70 AU and masses 
~ 0.002 — 0.07 Mq, within only 40 kyr. A comparison with 
observations shows that there are discs with similar proper- 
ties (see Section |4]). However, it is not possible to establish 
with certainty that fragmentation has occurred. The pres- 
ence of companions may be indicative of fragmentation, but 
these companions (i) may be too close to the central proto- 
star to be resolved as separate objects, (ii) may be too far 
away from the central protostar (e.g. within a few thousand 
AU) to be observed and/or associated with it, and (iii) may 
be below the sensitivity limit if they are not attended by 
an envelope/disc with enough mass to emit considerably in 
submm/mm wavelengths. 

Discs that do not fragment dissipate on a viscous 
timescale, i.e. on the order of 1 Myr. For example the non- 
fragmenting discs in our simulations (Fig. 3 and Table [l]) 
lose only a small fraction of their mass (as material is ac- 
creted onto the protostar) and increase in size as angular 
momentum is redistributed within the disc. Hence, if mas- 
sive (J; 0.2 Mq), extended {}i 100 AU) discs do not fragment 
then they should be readily and routinely observed, which 
is not the case. Therefore such discs either form and quickly 
fragment after becoming unstable due to mass loading from 
the surrounding envelope, or that their formation is sup- 
pressed (e.g. by magnetic fields). 

We conclude that current disc observations cannot ex- 
clude the mechanism of disc fragmentation as a means of 
producing low-mass hydrogen burning stars, brown dwarfs, 
and planetary mass objects. Interferometric observations of 
larger samples of Class objects with, e.g., ALMA may 
however provide a much stronger test of this scenario in the 
future. 



ACKNOWLEDGEMENTS 

Simulations were performed using the Cardiff HFC Cluster 
Merlin. Colour plots were produced using splash (Price 
2007). DS and APW acknowledge support by the STFC 
grant ST/HH001530/1. The work presented in this pa- 
per was stimulated by discussions held in the context of 
the Marie Curie Research Training Network Constellation 
(MRTN-CT2006-035890) . 



REFERENCES 

Andre, P., Ward-Thompson, D., & Barsony, M. 2000, Pro- 
tostars and Planets IV, 59 

Andre, P., et al. 2010, A&A, 518, L102 

Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C, & 
DuUemond, C. P. 2009, ApJ, 700, 1502 

Attwood, R. E., Goodwin, S. P., Stamatellos, D., & Whit- 
worth, A. P. 2009, A&A, 495, 201 

Balsara, D. S. 1995, Journal of Computational Physics, 121, 
357 

BeU, K. R., & Lin, D. N. C. 1994, ApJ, 427, 987 
Black, D. C, & Bodenheimer, P. 1975, ApJ, 199, 619 



10 Stamatellos et al. 



Boley, A. C. 2009, ApJ, 695, L53 

Boss, A. P., & Bodenheimer, P. 1979, ApJ, 234, 289 

Boss, A. P., & Myhill, E. A. 1992, ApJS, 83, 311 

Brinch, C, Crapsi, A., j0rgensen, J. K., Hogerheijde, 

M. R., & Hill, T. 2007, A&A, 475, 915 
Brinch, C, J0rgensen, J. K., & Hogerheijde, M. R. 2009, 

A&A, 502, 199 
Cai, K., Durisen, R. H., Michael, S., Boley, A. C, Mcji'a, 

A. C, Pickett, M. K., & D'Alessio, P. 2006, ApJ, 636, 

L149 

Clarke, C. J. 2009, MNRAS, 396, 1066 

Commercjon, B., Hennebelle, P., Audit, E., Chabrier, G., & 

Teyssier, R. 2010, A&A, 510, L3 
Cossins, P., Lodato, G., & Testi, L. 2010, MNRAS, 407, 

181 

Crutcher, R. M., Nutter, D. J., Ward-Thompson, D., & 

Kirk, J. M. 2004, ApJ, 600, 279 
Eisner, J. A., Hillenbrand, L. A., Carpenter, J. M., & Wolf, 

S. 2005, ApJ, 635, 396 
Eisner, J. A., & Carpenter, J. M. 2006, ApJ, 641, 1162 
Eisner, J. A., Plambcck, R. L., Carpenter, J. M., Corder, 

S. A., Qi, C, & Wilner, D. 2008, ApJ, 683, 304 
Enoch, M. L., Evans, N. J., Sargent, A. I., & Glenn, J. 

2009a, ApJ, 692, 973 
Enoch, M. L., Corder, S., Dunham, M. M., & Duchene, G. 

2009b, ApJ, 707, 103 
Evans, N. J., et al. 2009, ApJS, 181, 321 
Falgarone, E., Troland, T. H., Crutcher, R. M., & Paubert, 

G. 2008, A&A, 487, 247 
Forgan, D., Rice, K., Stamatellos, D., & Whitworth, A. 

2009, MNRAS, 394, 882 
Gammie, C. F. 2001, ApJ, 553, 174 

Goodman, A. A., Benson, P. J., Fuller, G. A., & Myers, 

P. C. 1993, ApJ, 406, 528 
Greaves, J. S., Richards, A. M. S., Rice, W. K. M., & 

Muxlow, T. W. B. 2008, MNRAS, 391, L74 

Hcnnebollo, P., & Ciardi, A. 2009, A&A, 506, L29 
Hennebelle, P., & Fromang, S. 2008, A&A, 477, 9 
Hennebelle, P., & Chabrier, G. 2008, ApJ, 684, 395 
Hennebelle, P., & Teyssier, R. 2008, A&A, 477, 25 
Rubber, D., Batty, C, McLeod, A., Whitworth, A., 2011 

A&A submitted 
Hubeny, I. 1990, ApJ, 351, 632 

Inutsuka, S.-I., Machida, M. N., & Matsumoto, T. 2010, 

ApJ, 718, L58 
Johnson, B. M., & Gammie, C. F. 2003, ApJ, 597, 131 
J0rgensen, J. K., et al. 2005, ApJ, 631 
J0rgensen, J. K., van Dishoeck, E. F., Visser, R., Bourke, 

T. L., Wilner, D. J., Lommen, D., Hogerheijde, M. R., & 

Myers, P. C. 2009, A&A, 507, 861 
Krasnopolsky, R., Li, Z.-Y., & Shang, H. 2010, ApJ, 716, 

1541 

Li, Y., Klessen, R. S., & Mac Low, M.-M. 2003, ApJ, 592, 
975 

Lin, D. N. C, & Pringle, J. E. 1987, MNRAS, 225, 607 
Lodato, G., & Rice, W. K. M. 2004, MNRAS, 351, 630 
Lodato, G., & Rice, W. K. M. 2005, MNRAS, 358, 1489 
Lyndcn-BcU, D., & Pringle, J. E. 1974, MNRAS, 168, 603 
Mayer, L., Lufkin, G., Quinn, T., & Wadsley, J. 2007, ApJ, 
661, L77 

Machida, M. N., & Matsumoto, T. 2010, arXiv: 1008.0920 



Machida, M. N., Inutsuka, S.-i., & Matsumoto, T. 2010, 

ApJ, 724, 1006 
Machida, M. N., Tomisaka, K., Matsumoto, T., & Inutsuka, 

S.-i. 2008, ApJ, 677, 327 
Machida, M. N., Matsumoto, T., Hanawa, T., & Tomisaka, 

K. 2005, MNRAS, 362, 382 
Masunaga, H., & Inutsuka, S.-I. 2000, ApJ, 531, 350 
Matzncr, C. D., & Levin, Y. 2005, ApJ, 628, 817 
Maury, A. J., Andre, Ph., Hennebelle, P., Motte, F., Sta- 
matellos, D., Bate, M., BcUoche, A., Duchene, G., Whit- 
worth, A. 2010, A&A, 512, A40 
Mellon, R. R., & Li, Z.-Y. 2009, ApJ, 698, 922 
Mellon, R. R., & Li, Z.-Y. 2008, ApJ, 681, 1356 
Meru, F., & Bate, M. R. 2010, arXiv:1011.5896 
Morris, J. P., & Monaghan, J. J. 1997, J.Comp.Phys., 136, 
41 

Nakano, T., Nishi, R., & Umebayashi, T. 2002, ApJ, 573, 
199 

Nelson, A. F. 2006, MNRAS, 373, 1039 

OfTncr, S. S. R., Klein, R. I., McKee, C. F., & Krumholz, 

M. R. 2009, ApJ, 703, 131 
Osterloh, M., & Beckwith, S. V. W. 1995, ApJ, 439, 288 
Padoan, P., & Nordlund, A. 2004, ApJ, 617, 559 
Price, D. J. 2007, Publications of the Astronomical Society 

of Australia, 24, 159 
Rodriguez, L. F., Loinard, L., D'Alessio, P., Wilner, D. J., 

& Ho, P. T. P. 2005, ApJ, 621, L133 
Spiegel, E. A. 1957, ApJ, 126, 202 

Stamatellos, D., & Whitworth, A. P. 2009a, MNRAS, 392, 
413 

Stamatellos, D., & Whitworth, A. P. 2009b, MNRAS, 400, 
1563 

Stamatellos, D., & Whitworth, A. P. 2008, A&A, 480, 879 
Stamatellos, D., Whitworth, A. P., Bisbas, T., & Goodwin, 

S. 2007a, A&A, 475, 37 
Stamatellos, D., Hubber, D. A., & Whitworth, A. P. 2007b, 

MNRAS, 382, L30 
Terebey, S., Shu, F. H., & Cassen, P. 1984, ApJ, 286, 529 
Thies, I., Kroupa, P., Goodwin, S. P., Stamatellos, D., & 

Whitworth, A. P. 2010, ApJ, 717, 577 
Tobin, J. J., Hartmann, L., & Loinard, L. 2010, ApJ, 722, 

L12 

Toomre, A. 1964, ApJ, 139, 1217 

Vorobyov, E. I., & Basu, S. 2010, ApJ, 714, L133 

Walch, S., Naab, T., Whitworth, A., Burkert, A., & 

Gritschneder, M. 2010, MNRAS, 402, 2253 
Whitworth, A. P., & Stamatellos, D. 2006, A&A, 458, 817 
Zhu, Z., Hartmann, L., & Gammie, C. 2009, ApJ, 694, 1045 



